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Abstract 

We introduce a new approach to constructing analytic solutions of the 
linear PDEs describing elastodynamics. This approach is illustrated for 
the case of a homogeneous isotropic half-plane body satisfying arbitrary 
initial conditions and Lamb's boundary conditions. A particular case of 
this problem, namely the case of homogeneous initial conditions, was first 
solved by Lamb using the Fourier-Laplace transform. The solution of the 
f*"^ ■ general problem can also be expressed in terms of the Fourier transform, 

but this representation involves transforms of unknown boundary values. 
This necessitates the formulation and solution of a cumbersome auxiliary 
problem, which expresses the unknown boundary values in terms of the 
Laplace transform of the given boundary data. The new approach, which 
is applicable to arbitrary initial and boundary conditions, bypasses the 
' above auxiliary problem and expresses the solutions directly in terms of 

the given initial and boundary data. 

C3 _ 

Keywords: elastodynamics, half space, initial boundary value problem, Lamb's prob- 
lem, global relation. 



1 Introduction 

The problem considered in this paper has a long and illustrious history, which 
begins with the classic works of Sir Horace Lamb in 1904 [TJ[31[3]. In [T|, Lamb 
treated four basic problems, the so called Lamb's problems, including both 
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the two and the three dimensional cases. These problems have been studied by 
several authors, see for example [4j [5] [6] [7] HI [9j QUI [TTJ [12] . For the axisymmetric 
point-load problem, Pekeris 14. obtained an exact expression for the surface 
motion in the case that the point load is given by a Heaviside function. This 
result was further extended by Schiel and Protazio [B] . For the three dimensional 
cases, Payton, using the Betti-Rayleigh reciprocal theorem [5], obtained the 
displacement and the surface motion. Gakenheimer, in his PhD thesis [7], made 
a systematic investigation of the travelling normal point load, by using the 
Cagniard-deHoop method. Bakkeraq et al. in [S] revisited this problem and 
presented a straightforward implementation of the Cagniard-deHoop method. 
Several related problems of contact mechanics and wave propagation can be 
found in the books of Miklowitz [2 , of Johnson [3] , and of Payton [13] . 

The extensive study of Lamb's problems is perhaps due to the fact that 
they arise in a variety of physical situations, as well as to their interesting 
mathematical structure. 

In this paper, we study two-dimensional problems under the usual assump- 
tion of plane strain. We consider arbitrary initial conditions and general stress 
boundary conditions, including normal line load, tangential line load, and mixed 
line load. We refer to the latter three stress conditions as Lamb's boundary con- 
ditions. 

Most studies of Lamb's problems are based on the Helmholtz decomposition 
and on the use of the Laplace transform in time. In particular, Lamb was the 
first to use the above approach, as well as a certain superposition which is in 
fact equivalent to the use of the Fourier transform. Miklowitz made use of the 
Cagniard-deHoop method [2], which itself is based on the Laplace transform. 
Payton made use of Green's functions and of the Laplace transform. 

Helmholtz decomposition has the advantage of decomposing the P-waves 
and S-waves. However, it has the disadvantage of introducing higher order 
derivatives to the boundary conditions. The use of Laplace transform in time, 
essentially restricts the problem to the case of homogeneous initial conditions. 

Our new approach is based on the unified approach for solving linear and 
integrable nonlinear PDEs introduced by one of the authors in 14 . A crucial 
role in this method is played by the so-called global relations, which are algebraic 
relations coupling appropriate transforms of the unknown boundary values with 
transforms of the given data. For linear PDEs this method uses three novel steps 
[15]- [44]: 1. Derive a representation for the solution in terms of an integral 
involving a contour in the complex Fourier plane. This representation is not yet 
effective, because in addition to transforms of the given initial and boundary 
conditions, it also contains transforms of unknown boundary values. 2. Analyse 
certain transformations in the complex Fourier plane which leave invariant the 
transforms of the unknown boundary values. 3. Eliminate the transforms of the 
unknown boundary values, by combining the results of the first two steps and 
by employing Cauchy's theorem (or more precisely Jordan's lemma). 

This paper is organised as follows. In section 2, we recall the governing 
equations of elastodynamics and derive the global relations. In section 3 we 
implement step 1 and in section 4 we implement steps 2 and 3. In section 5 we 
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use the general representations derived in section 4 in order to analyse Lamb's 
problems. These results are further discussed in section 6. 

For certain complicated boundary value problems it seems that it is not pos- 
sible to eliminate from the integral representation of the solution the transforms 
of the unknown boundary values. However, for some of these problems, by using 
the global relations, one can derive expressions for the Laplace transforms of the 
boundary values in terms of the given data [45 . A summary of this less effective 
approach is presented in the Appendix. It is interesting that for the particular 
case of zero initial conditions, the formulae presented in the Appendix reduce 
to the formulae first derived in the classical works of Lamb. 



2 Governing Equations and Global Relations 

The transient problem for two dimensional elastodynamics in the half plane 
with Lamb's boundary conditions is defined as follows: Let u = u(x,y,t), v = 
v(x, y, t), denote the displacements of a homogeneous isotropic half space body. 
The governing equations of motion without external body forces, are the Lamc- 
Navier equations: 

(A + 2p)u xx + (A + p)v xy + fiu yy - pu tt = 0, (2.1a) 

(J-v xx + (A + fi)u xy + (A + 2p)v yy - pv tt = 0, (2.1b) 

— oo < x < co, < y < oo, t > 0, 

where A, fi are the Lame constants and p denotes the density of the material, 
which without loss of generality is normalized to unity, i.e., p = 1. Let the initial 
conditions be denoted by 

u(x,y, 0) = u (x,y), u t (x, y, 0) = m(x, y), (2.2a) 

v(x,y,0) = v (x,y), v t (x, y, 0) = y), (2.2b) 

— oo < x < co, < y < oo. 

Let the stress boundary conditions be denoted by 

(u y + v x )(x,0,t) = gi(x,t), — oo < x < oo, t > 0, (2.3a) 

( v v + *\ 7T~ u x) i x t 0) *) — 92{x,t), — oo < x < oo, t > 0. (2.3b) 

V A -f- 2p / 

For a tangential line load, the functions g\ and g2 are given by 

9l (x,t) = a 6{x)X(t)/p, g 2 (x,t) = 0; (2.4) 
for a normal line load, 

9l {x,t)=0, g 2 (x,t) = a 5(x)Y(t)/(X + /i); (2.5) 
for a moving normal line load with a constant speed C , 

gi{x,t) = 0, g 2 (x,t) = a 6(x-Ct)/(X + fJ,). (2.6) 

Here <7o is a constant, 8(x) denotes the Dirac-<5 function, X(t) and Y(t) are 
functions which depend only on t. 
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Notations Hat, "A" , will denote the two-dimensional Fourier transform with 
respect to x and y, whereas tilde, "~" , will denote the Fourier transform with 
respect to x. In particular, 

/OO /"00 
dx dy e~ ikx - ay u(x,y,t), (2.7a) 
- oo JO 

/oo />oo 
dx / dy e- ikx - ilv v(x,y,t), (2.7b) 
-oo JO 

/oo />oo 
dx / dy e- ikx - ily Uj (x,y), (2.7c) 
-oo JO 

/oo />oo 
/ dy e- lkx - llv v 3 {x,y), (2.7d) 
-oo JO 

fc e R, Z e C", i > 0,j = 0,1, 

where C _ denotes the lower half complex Z-plane. 
Furthermore, 

/oo 
dx e- ikx u(x,0,t), (2.8a) 
- OO 
p OO 

v(k,t);= dx e- tkx v(x,0 7 t), fceR,i>0. (2.8b) 



gj{k,t):= dx e- lkx gj(x,t), k G R, i > 0, j = 1, 2. (2.9) 

J — oo 

We emphasise that since a; G M, the x— Fourier transform is well-dehned 
only for k £ R; on the other hand, since < y < oo, the Fourier transform 
is well-defined for I in the lower half complex /—plane. 

By applying the two-dimensional Fourier transform to equations (|2.ip and 
by using in the resulting equations the boundary conditions (|2.3p , we obtain the 
following equations: 

— (A + 2/j,)k 2 u + (A + fji)ik{ilv — v) + fj.\ — l 2 ii — ilu — (gi — ikv)\ = Uu, 

(2.10a) 

— fik 2 v + (A + fj,)ik(ilu — u) + (A + 2/i) — l 2 v — ilv — (c/2 — : ^— ikii) 

= v tt , k £ R, I £C. 

(2.10b) 



Introducing the notations 



P(k, I, t) = ku(k, I, t) + lv(k, I, t), (2.11a) 
Q(k, I, t) = lu(k, I, t) - kv(k, I, t), (2.11b) 
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(2.13b) 



equations (j2.10j) become 

P tt + (\ + 2(i)(k 2 + l 2 )P = F P , (2.12a) 

Q tt + Kk 2 + i 2 )Q = f q , feel, leC, (2.12b) 

where the functions Fp(k,l,t) and Fq(k,l,t) are defined as follows: 
F P (k, l,t)=- \pkgi(k, t) + iAfc 2 w(fc, t)] — + 2^)g 2 (k, t) 

(2.13a) 

+ 2ifiku(k, t)] - il 2 (X + 2fi)v(k, t), 
F Q (k,l,t) =[fc(A + 2/i)<? 2 (M) + i^k 2 u(k,t)} - l\p,gx(k,t) 
— 2ifxkv(k, t)] — i/il 2 u(k, t). 

Solving equations (|2. 12|) for {P, Q}, we obtain the following expressions: 

ku + lv=^—( e - lullt [ e luJlS Fp{k,l,s)ds-e Wlt f e- lu > lS F P (k,l,s)ds) 
-2iu\ V Jo Jo 1 

(2.14a) 

lu- kv = — !— (e- lLJ2t f e luJ2S F Q (k,l,s)ds-e luJ2t [ e~ lul2S F Q (k,l, s)ds] 
-2iu 2 V J J J 

+(^+Sk-'+(^-s; 



(2.14b) 

k G K, leC, 

where the dispersion relations u± and uj 2 are given by 

Luf = (A + 2/i)(fc 2 + ; 2 ), ul = fi(k 2 + l 2 ), fceK, leC, (2.15) 

and the known functions Po(k, I), P\{k, I), Qo(k, I), Qi(k, I) are given in terms of 
the initial conditions by 

Po = kiiQ + Ivq, P\ = kiii + Ivi, (2.16a) 
Qo = Iuq — kvo, Pi = lu\ — kvi, (2.16b) 



with uj, 6j, j = 0, 1, defined in equations (|2.7p . 
In what follows, we take 

oji = v / A + 2/i(fc 2 + / 2 )2, uj 2 = Jji(k 2 +l 2 )i. (2.17) 

The function (k 2 + l 2 )i has the branch points ±ik in the complex I— plane; we 
connect these two branch points by a branch cut and we fix a branch in the cut 
plane by the requirement that, 

(k 2 + l 2 )^ ~ l + o( 7 ), as l^oo. 
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Let uw)± ) j>Cj')± ) £/u), yw, j = 1,2, denote the following unknown func- 
tions: 



u ij)± (k,l,t)= I e ±l ^ s u(k,s)ds, (2.18a) 
Jo 

wW^M.t) = / e ±l ^ s {i(fc,s)ds, (2.18b) 
Jo 

U^(k,l,t) = —(e- iu ^u^ + {k,l,t) - e^u^-ikj,^), (2.18c) 

V^(fc,i,*) = ^-(e- lu ^v {j)+ {k,l,t)^e lu ^v {3] -{k,l,t)), (2.18d) 

fc G M,i G C,i > 0,j = 1,2. 

Similarly, let g^*, /W ± , j = 1,2, denote the following bown 

functions: 



g^ ± {k,l,t) = e ± ^ s ~g 1 (k lS )ds, (2.19a) 
Jo 

/« ± (fc,2,t)= / e^^Mjds, (2.19b) 
Jo 

G«>(M,*) - ^-(e-*-^ + ( fc>Zjt )_ e ^* fl O)-( fc>Z>t )) ) (2.19c) 

F^(k,l,t) = ^- (e-^'/ (i)+ (^ M) - e fajt / (3V M)) , (2.19d) 

fe G M.,1 G C,t > 0,j = 1,2. 
Using the above notations, equations (|2.14[) become 

ku + Iv =(\k 2 + l 2 (X + 2(i))V {1) + 2nklU {1) + N P , (2.20a) 
lu~kv = - [i(k 2 - l 2 )U (2) - 2fiklV {2) + N Q , (2.20b) 

where the known functions Np(k,l,t) and Nq{k,l,t) are defined as follows: 

N P (k, l,t)=- il(X + 2^)F {1) [k, I, t) - ifikG (1) (k, I, t) 

+( ^ (t , + i^o )e -~, + ( ^ (M) _^ )e ~,, 

(2.21a) 

N Q (k,l,t) =ik(X + 2fi)F^(k,l,t)-inlG^(k,l,t) 



+ (W*.o + |2^2).- + (i<WM-i2f2) 

(2.21b) 

We will refer to equations (|2.20p as the global relations. The global relations 
express the two-dimensional Fourier transforms of the solution (u, v) in terms 
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of the given initial and boundary data, as well as in terms of the transforms 
[/( 2 ) ; of the unknown boundary values. The important observa- 

tion is that equations (|2.20l) are valid for all values of I in the lower half complex 
I— plane. It turns out that using this fact it will be possible to eliminate the 
unknown transforms. 



3 An Integral Representation Involving the Un- 
known Transforms 

We first observe that for any fixed k £ M. and any fixed t, < t < T, T > 0, the 
functions j £ {1, 2}, are analytic in the complex /—plane. 

Indeed, these functions are all single-valued, thus it only remains to estab- 
lish the analyticity in the neighbourhood I = ±ik. The definition of i.e. 
equation (|2.18cj) . implies that this function possesses the following expansion 
for I near ±ik: 

oo 2n-2 r-t 

[/« (k, I, t) = iJ2 7Tj / (s - t) 2n ~ 2 u(k, s)ds. (3.1) 

^ (2n - 1)! Jo 

Similar expansions are also valid for U^\V^\V^\G^\G^ 2 \F (1 \F^ . 
Solving equations (12.20[) we find 



u = tt^I (Afc 3 + kl 2 (\ + 2fi))V w + 2ak 2 lU w 

fc+PL (3.2a) 
- [i(k 2 l - l 3 )U {2) - 2ukl 2 V {2) + kN P + IN Q \, 

* = TT~~~F? \ i^ 2 1 + l 3 (X + 2n))VW + 2ukl 2 U^ 

fc+PL (3.2b) 
+ /i(fc 3 - kl 2 )U {2) + 2ak 2 lV {2) + IN P - kN Q y 

We observe that the following important transformation is valid in the com- 
plex I— plane: 

I —I maps uji to — lji, Ll>2 to — u>2, and 
leaves U^'jV^', j = l,2 invariant. 

Employing the transformation I — > — I in equations p.2a[) and (|3.2b[) and then 
adding the resulting equations to (I3.2a|) and (|3.2bl) . we obtain the following 
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equations: 

u(k,l,t) +u(k,-l,t) = 

1 - T L 1 r 2 {2(\k 3 + kl 2 {\ + 2fjL))V {1 \k,l,t) - ^kl 2 V^(k,l,t) + (3.4a) 

k(N P (k,l,t) + N P (k,-l,t)) + l(N Q (k,l,t) - N Q (k,-l,t))}, 
v(k,l,t) +v(k,-l,t) = 

^L^^kl 2 U^(k,l,t) + 2 f i(k 3 -kl 2 )U^(k,l,t)+ (3 . 4b) 
l(N P (k, I, t) - N P (k, -I, t)) - k(N Q (k, I, t) + N Q (k, -I, t))}. 
Applying the inverse Fourier transform formula to these equations, we obtain 
u{x,y,t) = —l dk I dl e lkx+lly u(k,l,t) 



— oo J — oo 



i poo poo ikx-\-ily , 

— / dk dl— -\2(Xk 3 + kl 2 (X + 2fi))V w (k,l,t)-4fikl 2 V i2) (kJ,t) 

+ fc(iV P (fc, I, t) + 7V P (fc, -I, *)) + l(N Q (k, I, t) ~ N Q (k, -I, *))}, 

(3.5a) 

v(x,y,t) = — dk dl e lkx+ily v(k,l,t) 



4tt 2 _ 

i p oo />oo ikx+ily , 

= 4^i * 7 dZ fc 2 + I 2 | 4 ^ 2[/(1) ( fc - *) + 2 (^ fc3 - likl 2 )U {2) (fc, I, t) 
+ Z(7V P (fc, Z, t) - N P (k, -I, t)) - k(N Q (k, I, t) + N Q (k, -I, i))}, 

(3.5b) 

— oo < x < oo, < y < oo, t > 0. 

Denote by Hi(k,l,t) and H2(k,l,t) the following functions appearing in 
equations ()3.5|) : 



Hi(k,l,t) = ,, 1 ,. 7 (2(Afc 3 + Zd 2 (A + 2^))yW(fc,Z,t) - VdV 2 )(fc,M)l, 
fc 2 + Z 2 I J 

(3.6a) 

H 2 (k,l,t) = , „ 1 79 (^fcZ 2 ^^) (fc, t, t) + 2(^fc 3 - ukl 2 )UW (fc, Z, t) \. (3.6b) 
fc 2 + Z 2 L ' J 

We observe that for any fixed fc € R and fixed Z, < Z < T, T > 0, the above 
functions are analytic in the entire complex Z— plane. Indeed, equation (|3.ip and 
the analogous equation for U^ 2 \ imply that in the neighbourhood of I = ±ik, 
the following expansion is valid: 

H 2 (k,l,t) = 2i[ik [ u(fc,s)ds + o(w 2 ). (3.7) 



8 



Similarly, 



H 1 {k,l,t) = 2i\k I v(k,s)ds + o(u}f). (3.8) 



o 



The restriction y > 0, as well as the analyticity of the functions l,t) 
and H2(k,l,t), allow us to deform the contour of integration from the real axis 
to a contour 7^ in the upper half /— plane (the particular choice of jk will be 
determined in the next section): 



u(x,y,t) 



dk / dl 



jU2 1 72 

-00 Jj k " ~ t 

dk I dl 



dk I dl e 



|2(A/c 3 + fc/ 2 (A + 2n))V {1 \k, I, t) - 4/ifc/V (2) (fc, I, t)\+ 

h 2 +l 2 {k[N P (k,l,t) + N P (k,-l,t)]+ 

l[N Q (k,l,t)-N Q (k,-l,t)]], 

(3.9a) 

v(x,y,t) = 

+l \^kl 2 U {1) {k,l,t) + (fik 3 - fj,kl 2 )U^(k,l,t)\+ 

roo ikx+ily f (3.9b) 

dk / * -T^—w\l[Np(k,l,t) - N P (k,-l,t)]- 

k[N Q (k,l,t) + N Q (k, -Z,i)]}. 

4 The Elimination of the Transforms of the Un- 
known Boundary Values 

Let 
and 

•--'(if* 



The function Z12 has two branch points i/c^/^j^, which we connect by 



a horizontal branch cut; the function I21 has two branch points ±iky -^^jj, > 
which we connect by a vertical branch cut, see FigfT] We fix the branches by 
the requirements that, 
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a. The branch cut for U 



b. The branch cut for Z21 • 



Figure 1: Branch cuts for l%% and Z21 in complex I— plane. 



The following transformations are valid in the cut /-plane: 

/ — » I21 maps ui2 to — to to (4.3a) 

I — > Z12 maps ui\ to — u>2, [/« to V™ to V<- 2 \ (4.3b) 



Using in equations (|2.20a[) and (|2.20b|) the transformations I — > I12 and 
I — > —I respectively, and then combining the two resulting equations, we obtain 
the following equations: 



ku(k, Z12, t) + h 2 v(k, h 2 , t) = (Afc 2 + l 2 12 (X + 2/i)) V {2] [k, I, t) 

+ 2 t ikl 12 U < - 2 ' ) (fc, I, t) + N P {k, l 12 , t), 

-lu(k,-l,t) - kv(k,-l,t) = -fi(k 2 - l 2 )U {2 \k,l,t)+ 

2jiklV {2) {k,l,t) + N Q {k,-l,t). 



(4.4a) 
(4.4b) 



Similarly, using in equations (I2.20a|) and (12.20b[) the transformations I — > — I 
and I — > I21 respectively, and then combining the two resulting equations, we 
obtain the following equations: 



ku(k,-l,t) - lv(k,-l,t) = (Afc 2 +l 2 (X + 2fi))V {1) (k 1 l 1 t) 

- 2fiklU {1) (k, I, t) + N P (k, -I, t), 
kiu(k, hi,t) - kv(k, hi,t) = -/i(fc 2 - l 2 2i)U (1) (fc, I, t) 

- 2 M fcZ 2i (fc, I, t) + N Q (k, hi,t). 



(4.5a) 
(4.5b) 
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Let 

d(k,l) = Xk 2 + l 2 12 (X + 2fi), C 2 (k,l)=2fikl 12 , (4.6a) 

C 3 (fc, = 2/iM, C 4 (fc, = -fi(k 2 - I 2 ), (4.6b) 

L>i(fc, = Afc 2 + l 2 (X + 2/i), L> 2 (fc, = -2/iM, (4.6c) 

L> 3 (fc, = -2fxkl 2 i, D 4 (k, I) = -^{k 2 - l 2 2l ), (4.6d) 

A x (A;,/) = Ci(fc,l)Gi(M) - C 2 (k,l)C 3 (k,l), (4.6c) 

A 2 (fc, - Z?x(fc, J)I>4(A, ~ £> 2 (fc, 0^s(fc, 0- ( 4 - 6f ) 

Simplifying the expressions for Ai and A 2 we find 

Ai(fc, = n 2 (k 2 - I 2 ) 2 - Afi 2 k 2 ll 12 , (4.7a) 
A 2 (fc, = (Afc 2 + l 2 (X + 2n)) 2 - 4n 2 k 2 ll 2 i. (4.7b) 

Equations (|4.4[) imply 

V (2) (k, I, t) = 1 [C 4 (fc, '12, *) + iia«(fc, I12, t) - JV P (*. Zia, t)) 

Ai(fc,i) 

- C 2 (k,l)(-lu{k,-l,t) - kv(k, -l,t) - iV Q (-0)], 

(4.8a) 

Z ' *) = A n n [- g 3(fc, 0(fc«(fc, iia, t) + /i2«(fc, Z12, t) - iV P (fc, ii2, *)) 

Ai(fc,i) 

+ C 1 (k,l)(-lu(k,-l,t) - kv{k,-l,t) - N Q (k,-l,t))]; 

(4.8b) 

Equations (I4.5[) imply 

(fc, /, i) = - ,J DAk, l)(ku(k, -I, t) - lv(k, -I, t) - N P (k, -I, t)) 
A 2 (k,l) 

- D 2 (k, l)(hiu(k, hi,t) - kv(k, l 21 ,t) - N Q {1 21 ))], 

(4.9a) 

U(1) Z > *) = a H ,A ~ D ^ k ' 0(*«(*. -i, *) - l HK -/, t) - ^Vp(fc, - J, *)) 
A 2 (fc,Z) 

+ £>i(fc, l)(hiu(k, h u t) - fcfi(fc, iai, t) - N Q (k, hi,t))\. 

(4.9b) 

We fix the choice of the contour by requiring that every term in the RHS 
of (|4.8I) and (|4.9p does not have a pole or a branch point above this contour, 
see FigH 

Regarding the zeros of Aj, j = 1, 2, we note that they are of the form I — ak, 
for some constant a € C. For example, if A = 2fi, we find that the zeros of Ai 
are 

I « (-1.624 ±0.126i)fc, / « ±0.357ifc, i « ±1.056ifc, lw (1.624 ± 0.126i)fc; 
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cut /-plane 



71- 



Figure 2: 7^, the deformed path of integration. 




whereas the zeros of A2 are 

I w (-0.295 ± O.U2i)k 7 I m ±0.885ifc, Z w ±ift, Z w (0.295 ± 0.442i)fc. 



Substituting equations (J478J) and (|4.9|) in equations (|3.9p , and using Jordan's 
lemma in the complex /—plane above the contour 7^, it follows that u(fc, — Z, t), 
u(k,li2,t), u{k,l%x,t), v(k,—l,t), v(k,li2,t), v(k,foi,t), yield a zero contribu- 
tion. Hence equations ()3.9j) become the following equations: 

u(x,y,t) = 

/j; ikx+ily , 
(fc2 + ^ 2)Ai(fc;0 {2(Afc 3 + fcZ 2 (A + 2 At ))[- J D 4 (fc, 0^Vp(fc, -M) 

2 (fc, Z)AT Q (fc, Z 2 i, i)] - ^kl 2 [-Ci(k, l)N P (k, l 12 , t) + C 2 (k, l)N Q (k, -I, *)] } 

d/c / d * I2 + ,1 \k[N P (k,l,t) + N P (k,-l,t)] + 
x J -00 k + 1 I 

l[N Q (k,l,t)-N Q (k,-l,t)}}, 

(4.10a) 

v(x,y,t) = 

j; ikx+ily , 

(kz + F)A 2 (k { 2 ^ 2 [-P3(fc, Z)iV P (fc, -Z, t) - Z)7V Q (fc, Z 21 , *)] 
( M fc 3 - fikl 2 )[C 3 (k, l)N P (k, l 12 ,t) - d(fc, l)N Q (k, -I, t)]} 
^J^{z[iVp(fc, Z, i) - 7V P (fc, -Z, *)]- 

*[JV Q (M,i) + -'»*)]}> 

(4.10b) 

where {Cj,Z>j}f and {A.,} 2 are defined in (I4.6[) and (|4.7|) . and the known func- 
tions Np and JVq are defined in ()2.21j) . 



dk 



dk / dl 

— OO —OO 
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We summarize the above result in the following proposition: 

Proposition 4.1 Let (u,v) satisfy the Lame-Navier equations (|2.1[) in the half 
plane 

-co < x < oo, < y < oo, t > 0, 

with the initial conditions (|2.2[) and i/ie stress boundary conditions (|2.3p . As- 
sume i/ia£ i/ie given functions 

{u j {x,y),v j (x,y)} j=0 ,i, {9j(x,t)}j=i,2 

have sufficient smoothness and decay. A solution of the above initial-boundary 
value problem, which decays for large (x,y), is given by equations (|4.10[) , where: 

a) The known functions (Np,Nq) are defined in (|2.21[) in terms of the trans- 
forms of the initial and boundary data (see equations (|2.15l) - (|2.19p ). 

b) {Cj,Dj}j and {Aj}1 are given by equations (|4.6|) and (|4.7I) . 

c) The contours depicted in Fig\^ are deformations of the real axis and 
determined by the requirement that the zeros of Aj, j = 1, 2 and the branch 
points of I21 and Z12 are all below 7^. 



5 The Normal Line Load with the Homogeneous 
Initial Conditions 

Consider a normal line load suddenly applied to an isotropic elastic half plane 
body. In this case, 

u = ui = v = vi = 0; gi = 0, g 2 = a S(x)h(t)/(X + fi), (5.1) 

where ctq is a constant and h(t) is the Heaviside function defined by h(t) = 
0, t < 0; h(t) = 1, t > 0. 

In what follows, we compute the functions needed in equations (|4.10p . Equa- 
tions (1231) and (l2~Tg|) imply 

gi{k, t) = 0, g 2 (k, t) = -r^-h{t), (5.2a) 

A + /I 

g W ± (k,l,t) = 0, G u) (k,l,t) = 0, (5.2b) 
= -£$_J_( e ±*** - 1), (5.2c) 

A + fj, ±lU!j 

F^W) = -^(±-^¥ 1 ), J =1,2. (5.2d) 
Thus, the known functions Np and Nq defined in (|2.21l) . are given by 

N P (k,l,t) = -la -—^ (-J -ijfli), (5.3a) 

/ , , \ , A + 2a / 1 cos fcjot) \ , . 

N Q (k,l,t) = ko —^(— ^ii). 5.3b 
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The above equations imply 

N P (k,l,t) + N P (k,-l,t) = 0, (5.4a) 

N Q {k,l,t) + N Q (k,-l,t) = 2ka ±±^(± - ^^1), (5.4b) 

N P {k, l,t)- N P {k, -I, t) = -2l* ^ i- 2 - C ^ll) , (5.4c) 

N Q (k,l,t)-N Q (k,-l,t) = 0, (5.4d) 

N P (k,h 2 ,t) = - l -fN Q {k,l,t), (5.4e) 

N Q (k,l 2 i,t) = -jN P (k,l,t). (5.4f) 

Substituting the above relations and equations (|4.6p and (|4.7p into equations 
(|4.10p . we find the following result: 

Proposition 5.1 Let (u,v) satisfy the Lame-Navier equations (|2.1[) in the half 
plane 

— oo < x < oo, < y < oo, t > 0, 

wii/i homogeneous initial condition and the normal line load boundary condition 
(|5.ip . The solution of this initial-boundary value problem, which decays for large 
(x, y), is given by 

u{x,y,t) = 

A + 2/1 f°° f jkx+ily 

^°a + /iL 4 (fc 2 + P)W (fc 2 - z 2 ) 2 - 4 M 2 fc 2 z; 12 ] ( 5 - 5a ) 

|2(Afc 3 Z + M 3 (A + 2/i))(l - cos (wit)) + 4///cZ 2 Zi 2 (l - cos (w 2 *))}, 
v(x,y,t) = 

A + 2^ f 00 „ /" „ e ikx+ily 

— o-q— / dk dl 



A + m7-oo J lk (fc 2 + Z 2 )[(Afc 2 +Z 2 (A + 2 M )) 2 -4 M 2 fc 2 ZZ 21 ] 
|2/zfc 2 Z 2 [2^ 21 + (Afc 2 + Z 2 (A + 2 M ))] 1 ~ C ° S 2 (u;if) 

- 2 M 2 (fc 4 - /fc 2 Z 2 )(fc 2 - Z 2 - 2ZZ 12 ) 1 ~ CO& 2 (cJ2t) | 



2o~q — — — — / dk 
A + fi 



oo ^ e , Mv 2 l_ C0S ( {Jlt ) fc2 l-cosM| 
-oo fc 2 +Z 2 l W? " W| J ' 

(5.5b) 



where (cji,oj 2 ) are defined by equations (|2.17p and Zi 2 ,Z 2 i are defined by equa- 
tions (|4.2p , (|4.ip respectively. 



Similar results can be obtained for the other Lamb's problems. 
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6 Conclusions 



The main result of this paper is the derivation of equations (|4.10l) . These equa- 
tions express the displacements (u(x, y, t), v{x, y, t)) in terms of an integral along 
the real line and integrals along contours 7^ of the complex I— plane; these 
integrals involve transforms of the given initial and boundary data. Indeed, 
equations (|4.10p involve the functions Np(k,l,t) and NQ(k,l,t), which are de- 
fined in equations (|2.2ip in terms of the Fourier transforms Po(fc, I, t), Qo(k, I, t), 
Pi(k, I, t), Qi(k, I, t) of the initial data (see equations (|2.16I0 . as well as in terms 
of certain known transforms G^(k,l,t), G^(k,l,t), F^(k,l,t), F^(k,l,t) of 
the boundary data (see equations (12.19|) ). 

The starting point of the derivations of equations ()4.10|) is the derivation of 
the global relations (|2.20l) . These equations are the direct consequence of the 
application of the two-dimensional Fourier transform and of the substitution in 
the resulting equations of the given initial and boundary conditions. Equations 
(|2.20p involve the known functions Np and Nq, as well as the unknown functions 
U^ 2 \ (these functions involve certain transforms of the unknown 

boundary values u(x,0,t) and v(x,0,t), see equations (|2.8p and (|2.18[) ). The 
elimination of the above unknown functions is achieved by the following steps: 
1. By employing the transformation I — » —I, which leave the unknown functions 
(CA 1 ), EA 2 ) , , V^ 2 ') invariant, and by utilising the analyticity properties of 
these functions, we obtain the integral representations (|3.9[) . These represen- 
tations involve an integral along the real k — axis and an integral along the 
contour 7^ of the complex I— plane. 2. Using the transformations I — > Z21 and 
I -> l 12 which map to U^\V^ to and U^> to J/( 2 ),V« to re- 
spectively, we express the unknown functions {U^\ U^, V^) in terms of 
known functions as well as in terms of certain unknown functions which however 
are analytic in a certain domain of the complex I— plane, see equations (|4.8[) and 
(|4.9[) . 3. Using (|4.8[) and (|4.9p in equations p.9p and employing Jordan's lemma 
we obtain equations (|4. 10|) . 

The main advantages of the new approach are the following: 

1. The new method provides an analytic solution of Lamb's problem with 
arbitrary initial and boundary conditions. 

2. This solution is expressed in terms of the given initial and boundary data. 
The relevant representation is novel even for the particular case of homo- 
geneous initial conditions (this case has been analysed by several authors). 
An alternative approach using the Laplace transform is briefly discussed 
in the Appendix. By comparing equations (|4. 10[) and equation (|6.ip . the 
advantage of the new formulae becomes clear. Actually, taking into con- 
sideration that the initial-boundary value problems of the equations of 
elastodynamics are well posed for any finite t, the use of the Laplace 
transform, which requires t — > 00, is clearly inappropriate. 

3. The new method can be employed for the solution of several related initial- 
boundary value problems, including the problem of the isotropic half space 
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in the axisymmetric case. Furthermore, it can be extended to problems 
in three dimensions [46] . 



Appendix 



It is possible to analyse Lamb's problem by using only the transformations 
instead of using the transformations (|3.3p and the transformations (|4.3[) . How- 
ever, in this case one cannot eliminate directly all unknown boundary values. 
Instead, one can derive a complicated expression for the unknown boundary 
values in terms of the given initial and boundary data. (This approach is sim- 
ilar with the one used in |45) for solving Crighton's problem). Indeed, let the 



/-plane 




Figure 3: 7fc2, the path of integration. 

contour 7^2 be a simple curve in the lower half I— plane, determined by the 
requirement that it does not cross the branch cut associated with io\ and 102, 
see FigEJ Let K be the integral operator defined by 

™ )(M )- 1 r m, ' t)d ' 



2ir7,„ ((* 2 + P)'/2> 

for any function f(k,l,t) with appropriate smoothness and decay. Integrating 
the global relations (|2.20l) along 7^2 we find that the functions 



v(k,t) J ' "V""' V ft(M) 
satisfy a system of Volterra integral equations of the second kind: 

h(k, t) = (K[N})(k, t) * g(k, t) + (K[M]){k, t) * h{k, t) + (K[H])(k, t), (6.1) 

k e M, < t < T, 

where T > 0, * denotes the convolution operation with respect to t, 

( i^/JJle^ 2 * -i^/jl^j^ke lul2t \ 
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Mil, / +\ ( V**?^* ^^Jlkle^ \ 

\ VA+2u s/\+2n ' / 



eV" 2 * 2 v ^IfcZe^ 2 * 

2/J kj iuit A / 

and the known function H(k, I, t) is denned by 

H{k,l,t) = [i^-^-(ioj 1 P {k,l) + P 1 (k,l)), i— (iw 2 Q (fc,Z) + Q 1 (fc,0)J • 

The solution of the integral equations (|6.1[) yields the unknown transforms ap- 
pearing in (|2.20[) . 

Equations (|6.1[) can be solved in closed form by using the Laplace transform 
in t. Let "o" denote the Laplace transform with respect to t and let 

wi = (/ + (A + 2 M )/c 2 )5, W2 = (p 2 + M /c 2 )2, (6.2) 

Ag yy-.pij^, (6 . 3) 

V ^ + 2/1 W1W2 

The solution of (|6.ip is given by 



h = 



I VJlp 2 W2 ., (A+2M)(p 2 +2Mfc 2 )-2^/ A1 (A+2^)tt>itt> 2 

(p2 + 2^)2-4 v /^ M y II J 111 - 2 1 (p2 + 2 Atfc 2)2_4 v /3r Alfc 2 u)lU , 2 

■1 2 My/j^ : ^m'2-M(p 2 +2Mfc 2 ) yO+^p 2 -!!.! 



( P 2 +2 M fc 2 o- yT^fc \ / L f ^iPq+Pi 

^ U> 2 ][ 2tt J 7fc2 (A+2/i)(p-iwi) ;( fe2+/ 2)i 

_9,' / M p 2 +2pfc 2 II _L r W2Q0+Q1 rfi 

Z *Y A + 2^I tUl iJ| y \ 27T J 7fc2 H{p-iu 2 ) ;( fc 2 + ;2)5 

(6.4) 

The zeros of A coincide with the zeros of Rayleigh's function. When /i/A > 
0.906, the known transforms g(k,p), for each fixed fc, do noi have poles with 
positive real parts. 

In the particular case of the problem of homogeneous initial condition and 
the normal line load boundary condition, equation (|6.4p reduces to the classic 
Lamb's solution [T][2]. 
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